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Abstract 

We present a comprehensive study of two-boson exchange (TBE) corrections in parity- violating 
electron-proton elastic scattering. Within a hadronic framework, we compute contributions from 
box (and crossed box) diagrams in which the intermediate states are described by nucleons and 
A baryons. The A contribution is found to be much smaller than the nucleon one at backward 
angles (small e), but becomes dominant in the forward scattering limit (e — ► 1), where the nucleon 
contribution vanishes. The dependence of the corrections on the input hadronic form factors is 
small for Q 2 < 1 GeV 2 , but becomes significant at larger Q 2 . We compute the nucleon and A TBE 
corrections relevant for recent and planned parity- violating experiments, with the total corrections 
ranging from —1% for forward angles to 1 — 2% at backward kinematics. 



1 



I. INTRODUCTION 



Parity-violating electron-proton elastic scattering has become a standard tool with which 
to probe the strangeness content of the proton. Recent high-precision experiments at Jef- 



ferson Lab 



HQS 



4j and elsewhere 



8| have provided important constraints on the 
0, [k]] . Further improvements in the precision 



strange electric and magnetic form factors 
are expected to allow the measurement of the proton's weak char ge, Q m = 1 
where 6w is the weak mixing angle, to unprecedented accuracy fll. 12|. 



4 sin 2 9 



w ■ 



With the increasing precision comes the need to understand backgrounds to greater accu- 
racy than was called for in previous generations of experiments. In particular, higher-order 
radiative effects have received renewed attention recently, most notably those associated 



14 



16 



13, 



19 



20). 



with the exchange of two bosons (photons or Z-bosons) [13l . 
For point-like particles, the relevant loop diagrams are straightforward to compute and are 
included in the standard radiative corrections. However, incorporating the finite size of 
the nucleon leads to additional contributions, and can introduce further uncertainty in the 
calculations. 

In electromagnetic elastic scattering, despite being 0(a) suppressed, two-photon ex- 
change (TPE) was found to play an important role in resolving a large part of the dis- 
crepancy between the electric to magnetic proton form factor ratio measurements using the 
Rosenbluth and polarization transfer methods (see Ref. 2l| and references therein). One 
needs to carefully consider, therefore, to what extent the hadronic structure effects in two- 
boson exchange (TBE) may affect the analysis of parity-violating electron scattering. This 
is especially critical given that the extracted strange form factors appear to be rather small 
[9], as is the proton's weak charge Q w , which could further enhance the relative importance 
of TBE effects. 

In their seminal early work on electroweak radiative effects, Marciano & Sirlin [l^] com- 
puted the interference between the one-photon exchange and j-Z exchange amplitudes 
(which we denote by "7(^7)") at zero four-momentum transfer squared Q 2 , both at the 
quark level and at the nucleon level using dipole form factors. The corresponding contribu- 
tion from the interference between the single Z-boson and two-photon exchange amplitudes 
(denoted by "^(77)") vanishes at Q 2 = 0, but was computed within a generalized parton 
distribution formalism 16] at a scale Q 2 ~ several GeV 2 . 
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More recently, the 



including nucleon 



18 



BE corrections were computed at nonzero Q 2 in a hadronic basis, 



191 ] and A 20| intermediate states, with the structure dependence 



incorporated through hadronic form factors. For the nucleon intermediate states the model 



dependence was studied in Ref. 19[, and the individual TBE corrections to the proton and 
neutron terms in the parity-violating asymmetry computed. 

In this paper we perform a detailed analysis of TBE including both nucleon elastic and A 
intermediate states in the loop diagrams, and carefully examine their model dependence. We 



use the hadronic formalism developed in Ref. 22j], which allows a natural implementation of 
hadronic structure effects in radiative corrections at low Q 2 , where parity- violating electron 
scattering experiments are typically performed. For the A contribution we extend the two- 
photon exchange calculation of Kondratyuk et al. to the weak sector, and constrain the 
axial- vector form factors by data from neutrino scattering. 

In Sec. [Til we review the basic formalism of parity- violating electron scattering and sum- 
marize the Born level amplitudes and cross sections. The two-boson exchange corrections 
are described in Sec. IHIl where we outline the box diagram calculations with nucleon and A 
intermediate states. Our main results are presented in Sec. IIVI We compute the corrections 
from TBE to the parity-violating asymmetry, and discuss the consequences for the extrac- 
tion of the proton's strange form factors and weak axial charge. Finally, we summarize our 
findings in Sec. |V] and identify possible future developments of this work. 

II. BORN APPROXIMATION 

For elastic scattering of an electron e~ from a nucleon N we define the initial e~ and iV 
momenta as pi and P2, and final e~ and N momenta as ^3 and p^, respectively, e~(pi) + 
N(j>2) e ~{Pz) + N(p4). The four-momentum transferred from the electron to the nucleon 
is given by q = p^ — P2 = p\ — p$, with Q 2 = —q 2 > 0. In the Born approximation, the 
amplitudes for the electromagnetic and weak neutral currents are given by: 

~2 J-y "IP- ' 



g 2 1 G 



Mz ~ "(4cosM 2 W^ 2 3 ' ^ 3 " ' (2) 

where e is the electric charge, g = e/ sm9w is the weak coupling constant, Mz is the Z 
boson mass, and Gp = Tta./{y2M% sin 2 9 W cos 2 9 W ) is the Fermi constant, with a = e 2 /47r 



the fine structure constant. At tree level the weak mixing angle is related to the weak boson 
masses by sin 2 Q\y = 1 — M^/M§, where Mw is the W boson mass (in our numerical results 
below we use the renormalized value sin 2 9w = 0.2312 24j). The matrix elements of the 
electromagnetic and weak leptonic currents are given by 

# = uMYuM , (3) 

fz = u e (p 3 ) (g e v r + giris) « e (pi) , (4) 

where the latter is given by a sum of vector and axial-vector terms. We use the convention 
in which the vector and axial- vector couplings of the electron to the Z boson are given by 

g e v = -(l-4sin 2 M , g\ = +1 . (5) 

The matrix elements of the electromagnetic (weak) hadronic currents can be written as 

J 7(Z) = ^Ar(P4) U N (p 2 ) , (6) 

where the current operators are parameterized by the electromagnetic and weak form factors: 

n; = r f? n (q 2 ) + ^fi n (q 2 ) , (7) 

Y% = 7 - F^(g 2 ) + ^ F™ (Q 2 ) + 7 "75 Gr (Q 2 ) , (8) 

with M the nucleon mass. Here F\ and F 2 are the Dirac and Pauli form factors, and Ga 
the axial form factor of the nucleon (N = p,n), for either the electromagnetic (7) or weak 
(Z) current. Usually one takes linear combinations of the Dirac and Pauli form factors to 
define the Sachs electric and magnetic form factors as 

G E (Q 2 ) = F 1 (Q 2 )-rF 2 (g 2 ) , (9) 
G M (Q 2 ) = F 1 (Q 2 ) + F 2 (g 2 ) , (10) 

where r = Q 2 /4M 2 . 

The differential cross section is given by the square of the sum of the 7 and Z Born 
amplitudes, 

— =(^-Y\M\ 2 (11) 
where the squared amplitude can be written as 

|7W| 2 = \M 7 + M z \ 2 = \M^\ 2 + 23? (M;M Z ) + \M Z \ 2 . (12) 
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The purely weak contribution \Aiz\ 2 is small compared with the other terms and can be 
neglected. By polarizing the incident electron and measuring the difference between right- 
and left-handed electrons scattering from unpolarized protons, the parity-violating (PV) 
asymmetry can be defined in terms of the differential cross sections as 

A P y = , (13) 

where u r(l) is the cross section for a right- (left-) hand polarized electron. The purely 
electromagnetic contribution cancels in the numerator, so that the asymmetry is sensitive to 
the parity- violating part of 23? [Ai* Aiz) , involving the interference of .M 7 with the product 
of vector and axial-vector currents in Aiz (the vector-vector and axial-axial parts of Aiz 
cancel in the asymmetry). The denominator is dominated by the electromagnetic term, 
|.M 7 | 2 . 

More explicitly, the PV asymmetry can be written in terms of the electroweak form 
factors as 

where e and e' are kinematical parameters, 



iZN 



e' 1 = 1 + 2(1 +r) tan 2 - , (15) 



e' = y/ T (l + T)(l-e*) , (16) 

with 9 the electron scattering angle in the target rest frame. 

For a proton target the weak electric (magnetic) vector form factor G^ M ^ can be related 
by isospin symmetry to the electromagnetic form factors of the proton and neutron by 

G Z E{ M ) = (1 - 4 sin 2 Ow)G^e(,m) ~ ^E{M) ~ ^E(M) j 

where G% M -> are the contributions from strange quarks. The small factor (1 — 4sin 2 6 l jy) 
suppresses the overall contribution from the proton electromagnetic form factors, thereby 
promoting the neutron form factors to play a greater role. The weak axial-vector form factor 
of the proton is given by G 7 ^ = —G P A + G% where G S A is the strange quark contribution. 

Measurement of the PV asymmetry A-py as a function of the scattering angle 9 allows 
one to extract combinations of the strange form factors, given knowledge of the proton 
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and neutron electromagnetic form factors. Reliable extractions of the form factors require 
precise knowledge of the radiative corrections to the PV scattering associated with higher 
order electroweak processes. This is especially critical given that the extracted strange form 
factors appear to be rather small numerically. In the next section we discuss a subset of the 
radiative corrections, namely those arising from two-boson exchange. 



III. TWO-BOSON EXCHANGE CORRECTIONS 



Beyond the Born approximation, the PV asymmetry receives corrections from higher 
order radiative effects, such as vertex corrections, wave function renormalization, vacuum 
polarization, and inelastic bremsstrahlung, which are well known and included in standard 
data analyses. Less well determined are radiative corrections arising from the interference 
of Born and TBE diagrams, both electromagnetic (77) and electroweak ilZ). For purely 
electromagnetic scattering, the TPE corrections these have been shown 2l|, |22J to display 
strong angular dependence, which significantly affects extractions of the G^/G]^ ratio by 
Rosenbluth separation [25]. 



There are several ways in which the PV asymmetry can be represented in the presence 
of higher-order radiative corrections. The approach pioneered by Marciano & Sirlin [jjj] 
parameterizes the electroweak radiative effects in terms of parameters p and k, such that 
the weak charge of the proton in the presence of higher order corrections becomes 



1-4 sin 2 9 W 



p(l - 4k sin 2 6 



(18) 



In this case the asymmetry can be written as a sum of proton vector, strange vector, and 
axial-vector contributions, 

G F Q 2 



where 



.4 



PV 



AV2- 



TXOl 



(A v + A S + A a) 



Av = 9aP 



'1 - An sin 2 6 W ) 



Crcd 



1 

C"red 



A s = -g% P — (eGfG E + tG^G s m ) , 



A a — 9v £ ' G^G^i 



Orcd 



(19) 



(20a) 
(20b) 
(20c) 
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with cr rc d = eiG 1 ^) 2 + t(G]^) 2 the reduced unpolarized proton cross section. 

An alternative parameterization is in terms of isoscalar and isovector weak radiative 
corrections for the vector form factors, and a similar set of corrections for the axial-vector 
form factors. In this case the vector part of the PV asymmetry is written 



A 



v 



91 



1 - 4 sin 2 9 W ){\ + R P V ) (sGfG^ + tG^GZ) (1 + R 



0"red 



(21) 



where the proton and neutron radiative corrections are given, to first order in p — 1 and 
k - 1, by 

, 4 sin 2 9 W 



r?" 

tiy 



P-1- 
P-I. 



IK 



4 sin 2 W 



(22a) 
(22b) 



The strange part of the asymmetry, 

„ 1 



As 



-9a 



Orcd 



(sGJG% + rGZG s M ) (1 + R 



(23) 



receives an isoscalar radiative correction, given by 



R 



(0) 



p-1 . 



(24) 



For the axial asymmetry A A , the form factor 5? implicitly contains higher order radiative 
corrections for the proton axial current, as well as the hadronic anapole contributions |9l. 1261] . 



At tree leve 
In Refs. 



, and in the absence of the anapole term, G z ^ — > G^ p - 



18 



19 



201 ] the contributions to p and k from the interference of the Born and 



TBE (box and cross-box) diagrams were computed, denoted by Ap and An, respectively. 
The correction to the PV cross section arising from the the 77 and 7Z TBE contributions 
can be obtained from Eq. ffl2l by the replacements 



M- 



My + M, 



l 7 J"l-y I S"l-yy , 

M z -> M z + M yZ + M Zj 



(25a) 
(25b) 



where the two-photon and 7Z exchange amplitudes .A/f 77 , Aiz-y and .M 7 z are given explicitly 
below. The relative corrections from the ^(77), 7(7^), and 7(77) interference terms can 
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be identified as 



^(ty) 



2$t{M* z M 



id 
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7(77) 



2^{M Z M 1 ) ' 

2^{m* 1 m iZ + m;m Zi ) 

2^{M*M z ) 

2 k (A4;a4 77 ) 



(26a) 
(26b) 
(26c) 



The correction to the Born level PV asymmetry A PY can then be represented as 
A PV = (l+8)A% EE 



1 + £z(77) + &y(Z 7 ) A -.0 

1 + t) 7(ry) J 

where Apy is the full asymmetry, including TBE corrections, and A PV is given in 
Since the electromagnetic TPE correction <5 7 ( 77 ) is typically only a few percent 
the full correction 8 can be written approximately as 

8 ~ ^Z(77) + ^7(^7) — ^7(77) • 
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22 



(27) 



(28) 



In the model discussed here, the amplitudes .M 77 , M.^z and M. z ^ contain contributions 
from both nucleon elastic and A(1232) isobar intermediate states, which we discuss next. 

A. Nucleon Intermediate States 



For completeness, here we review the basic elements of the TBE exchange calculation with 



nucleon intermediate states. A more complete account can be found in Refs. 



3, 



21 



For electromagnetic scattering, the total 27 exchange amp 
diagrams with a nucleon intermediate state has the form 

d A k 



itude for the box and crossed-box 



2l| 



M 



77V7 



(2k)* 

x u N (p 4 ) F»(q - k) S F (p 2 + k, M) V v (k) u N (p 2 ) A F (k, A) A F (k - q, X) , (29) 



where m e is the electron mass, and the fermion (electron) and gauge boson (photon) prop 
agators are given by 

i (k + m) 



iSp(k, m) 
iA F (k,X) 



k 2 — m 2 + ie 
—i 

k 2 - X 2 + ie ' 



(30) 
(31) 



respectively, with A introduced as an infinitesimal photon mass to regulate the infra-red 
divergences. 

The calculation of the 7-Z interference amplitude proceeds along similar lines to that of 
the 27 amplitudes above, with the appropriate replacements of the photon propagator by 
the Z boson propagator, and the ~/NN vertex function by in Eq. (jSJ), 

eV f d A k 
M ^ NZ ~ (4cosM 2 iN I 

x u e (p 3 ) {gvln + 9a1^15)Sf(Pi ~ k, m)^ v + -f u S F (p 3 + k, m)(#y7^ + g e A ltil$) « e (Pi) 

x u N ( P4 ) r" z ( q - k) s F ( P2 + k, m) r;(k) u N ( P2 ) A F (k, \) A F (k - q , m z ) . (32) 

A similar expression holds for the conjugate amplitude M.zn^- 

For the electromagnetic nucleon form factors we use the global fit to the proton electric 
and magnetic form factors from Arrington et al. 25l, and for the neutron form factors 

n 

from Bosted [27J. For technical reasons, we parameterize the form factors by a sum of three 
monopoles. To examine the model dependence of the calculation, we also consider a dipole 
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22] 



shape for the proton form factors, with a dipole mass of Ajv(y) = 0.84 GeV 

The weak ZNN form factors are less well determined. Using the conservation of the 
vector current (CVC), the weak vector form factors can be directly related to the jNN 
form factors. For the axial-vector form factor, on the other hand, we use an empirical dipole 
fit, G A {Q 2 ) = GU(0)/(1 + Q 2 /A 2 N(A) ) 2 , where G A (0) = 1.267 is the axial vector charge, 
with the mass parameter An( A ) = 1 GeV. Varying An(a) by 20% does not affect the results 
significantly. Since the main purpose of the PV experiments is to extract strange quark 
contributions to form factors by comparing the measured asymmetry with the predicted 
zero-strangeness asymmetry, in all our numerical simulations we set the strange form factors 
to zero, F( 2 = = G S A . 

In Fig. [1] we show the various contributions to the two-boson exchange correction 5n 
'or several values of Q 2 (Q 2 = 0.01, 0.1, 1 and 5 GeV 2 ). The infrared 



1 of 


e : 


m 


29 



divergences 

Mo & Tsai [28(. It should be noted, however, that, in contrast to the 27 box diagrams, the 
infrared contributions for the 7Z box diagrams are significantly different using the procedure 
of Ref. 



291 ] . At small Q 2 values (Q 2 ^$ 0.1 GeV 2 ) the 7(77) and ^(77) contributions are 
very similar, and considerably smaller in magnitude than the 7(^7) component. Since the 
'y—Z interference and the purely electromagnetic contributions enter in the numerator and 



9 



0.04 



0.03 



O) 0.02 
to 



WD 



0.01 



.... y(yy) 

— Z(yy) 

— y(^y) 



Q~ = 0.01 GeV 2 



0.2 0.4 0.6 O.J 

e 






FIG. 1: TBE corrections 5n(£, Q 2 ) with nucleoli intermediate states, for the 7(77) (dotted), ^(77) 
(dashed) and 7(^7) (solid) contributions at Q 2 = 0.01, 0.1, 1 and 5 GeV 2 . The correction is 



defined relative to that of Mo & Tsai 
sign in the asymmetry, Eq. 



281 ] . Note that the 7(77) correction enters with the opposite 



denominator of the PV asymmetry, respectively, the 7(77) and ^(77) will partially cancel 
in their effect on Apy, which will be determined mostly by the 7(^7) component. At larger 
Q 2 ( > 1 GeV 2 ) the 7(^7) component decreases in magnitude, while the 7(77) ^(77) pieces 
become large and more negative 3, 
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22] . 



The dependence of the total correction Sn on the input form factors is illustrated in 
Fig. [21 The difference between the results using the empirical form factors and the dipole 
approximation is very small for all values of e, and only becomes appreciable at large Q 2 
(Q 2 ^ 1 GeV 2 ), consistent with the findings of our earlier analysis 21] . Interestingly, the 
correction at Q 2 = 0.01 GeV 2 is relatively flat over the range 0.1 < e < 0.8, before dropping 
rapidly as e — > 1. At large Q 2 the total TBE correction becomes more strongly e dependent, 
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0.04 
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0-02 
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8 

FIG. 2: Model dependence of the total TBE corrections 5n(£, Q 2 ) with nucleon intermediate states 
for Q 2 = 0.01, 0.1, 1 and 5 GeV 2 , using the empirical form factors as described in the text (solid), 
and using a dipole approximation (dashed). 

decreasing in magnitude at forward scattering angles but increasing at backward angles 
(e - 0). 




B. A Intermediate States 

In evaluating the contribution to the TBE amplitude from the excitation of the A(1232)- 
isobar, we use the formalism outlined in Ref. [23] for the "fNA interaction, and extend this 
to the weak sector with the introduction of axial ZNA couplings. The jNA vertex is given 



by [23|, [30| 



+ 92(Q 2 ) \f<t ~ <Tp ■ q] + 9j jP- [q 2 (p"7° - 9 m t) + <f (ft - 7> ■ ?)] ) 7s , 

(33) 

where p and q are the incoming A and photon momenta, with corresponding Lorentz indices 



a and /x, respectively. The overall factor \/2/3 arises from the N — > A isospin transition 
operator. Electromagnetic gauge invariance implies that q/jFt^A^iyip, q) = 0. The coupling 
constants s = QiiQ 2 = 0) for i = 1,2,3 can be related to the magnetic, electric and 
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Coulomb components of the ^NA vertex by g\ = <?m, gE = gi — <7i, = 93- The vertex 
with an outgoing A can be obtained from the relation 

F 7V->a(p> Q) = 7o [Fyl^ N (p, q)] f 7o , (34) 

where p is the outgoing A momentum and q the incoming photon momentum. 

The amplitude for the box and crossed-box diagrams with a A intermediate state can 
then be written as 

/d A k r 1 

■^yju e (p 3 ) i^Spipi - k, m e )^ v + j u S F (p 3 + k, m e )7^ w e (pi) 

x u N (p 4 ) T^_^ N (p 2 + k,q-k) S F (p 2 + k, M) V^{p 2 + k) r^_ A (p 2 + fc, fc) itjvfe) 
x A F (k,0) A F {k-q,0) , (35) 

where the projection operator 

Va^P) = 9aP ~ ^7a7/3 - T^irflaPp + VcTlfif) (36) 

ensures that only spin-3/2 components are present. Suppression of the unphysical spin- 1/2 
contributions also leads to the condition on the vertex p a T^^ N (p, q) = 0. Note that in 
Eq. (1351) a finite photon mass is not needed in the photon propagators, since, in contrast to 
Eq. (I29p . the result here is infra-red finite. 

For simplicity, we assume a dipole shape for the three 'yNA transition form factors, 
9i (Q 2 ) = g t Fft(Q 2 ) for i = 1,2,3, where Fft(Q 2 ) = (1 + Q 2 / A 2 A(y) )- 2 , with a dipole mass 
Aa(v) = 0-84 GeV for each. For the electric and magnetic couplings we use the values oi = 7 



and g 2 = 9 



obtained from a K-matrix analysis of pion photoproduction data 



30j. A 



more realistic ttN coupled channel quasi-potential study [3jJ gives similar values, g\ = 6.3 
and g 2 = 9.7. For the g% coupling, an estimate from the — > A E2/M1 transition strength 
yields #3 = 5.8. To test the sensitivity of the TBE corrections to the value of g^, we consider 
a range of couplings, as discussed below. Note that the interference contributions between 
the e/i, g 2 and g% terms cancel in the TBE amplitude because of the odd and even character 
of these vertices in the loop variable k. 

For the ZNA vertex both vector and axial-vector contributions enter. For the vector 



12 



transitions, CVC requires the same form for the ZNA vertex as for the 7 A" A, 



ZA— ►iV 



(p,q) 



2M| V 3 



{g v AQ 2 ) W°i i ~ P^?i - Yi a p ■ ? + 7V 



+ g%(Q 2 ) [pV - 9 m P ■ q] 



M A 



[q 2 (pV - gTj) - <f (qy - 7> ■ ?)] } 75 



(37) 



where again the factor \/2/3 is associated with the A" — > A weak isospin transition. Using 
CVC and isospin symmetry, the vector ZNA form factors can be related to the 7iVA form 
factors by 

^(Q 2 )=2(l-2sin 2 M^(Q 2 ), (38) 

where the Q 2 dependence of the electromagnetic 7iVA form factor is parameterized as above. 

For the axial-vector vertex, nonconservation of the axial current implies the existence of 
an addition form factor. However, one can use the partially conserved axial current (PCAC) 
hypothesis to relate two of the form factors, leaving a similar expression to that in Eq. (13?]) . 



ZA^N 



2M| 



{gt{Q 2 )[g^4 

9i(Q 2 



^ lr f L 4 - Yi a p • q + 7 /i <? Q i 



+ 9?(Q 2 ) [P"q a - 9 m P ■ q] + [q 2 (p M 7° - 9^) ~ <f W - TP ■?)]}• 



(39) 



Note that here the weak isospin transition factor has been absorbed into the definition of 



the couplings 



The axial form factors are less well determined, but some constraints 



have been extracted from analysis of v scattering data. In a recent analysis, Lalakulich & 



Paschos [32j parameterized the vN — > p,A cross sections from bubble chamber experiments 
at low Q 2 in terms of phenomenological form factors. The available data can be described 
by the form factors gf{Q 2 ) = 0, g£(Q 2 ) = {Ml/2M 2 )C£{Q 2 ) = (Q 2 /AM 2 ) g£{Q 2 ), where 



C5 is given in Appendix A, with C§(Q = 0) = 1.2 |32|. For the Q dependence we again 



take a dipole form, with a cut-off mass of Am a.) = 1.0 GeV. 

As for the electromagnetic case, the vertex with an outgoing A can be obtained from the 
relation 



7o 



v TZ\n{pa) 



t 



7o 



(40) 



where p is the outgoing A momentum and q the incoming Z-boson momentum. The ZNA 
amplitude for the box and crossed-box diagrams with a A intermediate state can then be 
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FIG. 3: TBE corrections 5 a (£, Q 2 ) with A(1232) intermediate states, for the 7(77) (dotted), ^(77) 
(dashed) and 7(^7) (solid) contributions at Q 2 = 0.01, 0.1, 1 and 5 GeV 2 . 



written 
M^az 



eV 



(4 cos 9w) 



d 4 k 



x u e (j> 3 ) {gvln + 9a1h1s) S f{Pi - k, m)^ v + ^ u S F {p 3 + k, m)(g-y7 jU + g e A ^ 5 ) u e {p x ) 

3/2 

»/3 



x «Ar(p 4 )r^ Ar (p 2 + k, q - k)S F (p 2 + k, M) V% 2 (p 2 + k) 



x t ^a(P2 + k, k)u N {p 2 ) A F (k, 0) A F (k - q, M z ) 



(41) 



where T Z \_^ N is the sum of the vector (137|) and axial- vector (139|) vertices. The corresponding 
amplitude -M 7 az can be derived in a similar manner. 

In Fig. [3]we plot the individual TBE contributions to <5a from processes with intermediate 
A(1232) states as a function of e for a range of Q 2 values between 0.01 and 5 GeV 2 . Several 
interesting features can be noted. Firstly, the magnitude and shape of the A corrections 
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are very different to the nucleoli corrections in Fig. [TJ At low Q 2 0.1 GeV 2 ) the two- 
photon interference with either the Born 7 or Z exchange is almost negligible, increasing 
somewhat at larger Q 2 . The j(Zj) contribution is also relatively small at low e, and none 
of the corrections exceed ~ 1% in magnitude for e < 0.8 and Q 2 < 1 GeV 2 , and ~ 2% for 
Q 2 < 5 GeV 2 . 

At larger e, however, the j(Zj) correction increases rapidly, becoming even bigger than 
the nucleon correction, and in fact appears to diverge as e — > 1. The increase of the one- loop 
contributions to the asymmetries may be related to the growth of the invariant center of mass 
energy for fixed Q 2 as s — > 1. Since the A intermediate state amplitudes .M 7 A7 and M^az 
have numerators which have higher powers of loop momenta than the corresponding nucleon 
amplitudes .M 7 tv 7 and A4^nz, one expects that the A contributions should grow faster with 
invariant energy than the nucleon. It is also interesting to observe the cusp behavior of the 
7(77) and ^(77) corrections at Q 2 = 0.1 GeV 2 around e = 0.6, the kinematics of which 
corresponds to the threshold point of the e-A channel. 

The combined TBE correction from A intermediate states is shown in Fig.H](a), for various 
input form factors. In general the behavior of the total correction 5a is quite dramatic at 
high e, with the magnitude increasing as e — > 1. The total correction for Q 2 < 0.1 GeV 2 
is positive for most e values, but changes sign to become negative at larger Q 2 . As for the 
nucleon case, the dependence on the input form factors is relatively weak for all Q 2 < 1 GeV 2 , 
whether one uses empirical form factors for the vector ^NN or ZNN vertices or a dipole 
approximation for all the form factors. Similarly, the dependence on the dipole cut-off 
masses Aa(v,a) for the jNA and ZNA vertices is small for the same Q 2 range, Fig. IU(b). 
The sensitivity to the input form factors becomes more appreciable at larger Q 2 , however, 
as the Q 2 = 5 GeV 2 results demonstrate. One should caution, though, that at momentum 
transfers of Q 2 ~ 5 GeV 2 or higher the reliability of a purely hadronic resonance description 
of the TBE process is more questionable. 

Finally, the dependence of 5a on the Coulomb coupling constant #3 is illustrated in 
Fig. El where the total correction at Q 2 = 0.01 and 1 GeV 2 is shown for = —2 23|, 



and 5.8 



30] . The results with g s = —2 and are almost indistinguishable, while using the 



preferred coupling g 3 = 5.8 gives slightly smaller contributions for most e. One can conclude, 
therefore, that the uncertainty in the Coulomb coupling should not affect the overall results 
or conclusions. 
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FIG. 4: Total TBE correction S A (e,Q 2 ) with A(1232) intermediate states for Q 2 = 0.01, 0.1, 1 
and 5 GeV 2 . (a) Comparison between using empirical nucleon form factors (solid) and a dipole 
approximation (dashed), (b) Dependence on the iV — > A transition form factors, using the 
standard cut-offs A A (y) = 0.84 GeV, A A (yi) = 1-0 GeV as described in the text (solid), and the 
modified cut-offs A A(y) = 0.68 GeV, A A(j4) = 0.8 GeV (dashed). 
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FIG. 5: Total TBE correction 5 A (e,Q 2 ) with A(1232) intermediate states for Q 2 = 0.01 and 
1 GeV 2 , with different Coulomb couplings g% = —2 (solid), (dotted) and 5.8 (dashed). 
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FIG. 6: TBE corrections 6 N (e, Q 2 ) for the nucleoli (dashed) and S A (e, Q 2 ) for the A(1232) (dotted) 
intermediate states, and the sum (solid), for Q 2 = 0.01, 0.1, 1 and 5 GeV 2 . 



IV. EFFECTS ON OBSERVABLES 



A comparison of the total TBE corrections with nucleon and A(1232) intermediate states, 
together with their sum, is presented in Fig. [6]for Q 2 = 0.01, 0.1, 1 and 5 GeV 2 . As observed 
in the previous section, at small e (e < 0.6) the TBE correction at Q 2 < 1 GeV 2 is dominated 
by the nucleon elastic contribution. At larger e the A plays an increasingly important role, 
and generally exceeds the nucleon piece at s > 0.9. At higher Q 2 , the magnitude of the 
A contribution is larger than that of the nucleon for most e values, although as remarked 
above, the reliability of a purely resonant description of TBE is less clear at momentum 
transfers above Q 2 ~ 5 GeV 2 . 

The Q 2 dependence is more clearly illustrated in Fig. [7J where we show the nucleon and 



17 



0.01 



0.1 1 
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FIG. 7: Total TBE corrections 5n (upper three curves) and <5a (lower three curves) versus Q 2 for 
fixed e values, e = 0.1 (dashed), 0.5 (dotted) and 0.9 (solid). 

A corrections for fixed e = 0.1, 0.5 and 0.9. At low Q 2 the nucleon correction 5n increases 
as Q 2 — * 0, but flattens out somewhat for larger Q 2 . The A correction <5a, in contrast, is 
almost Q 2 independent for Q 2 < 1 GeV 2 , except at very high e, but rapidly becomes large 
and negative at higher Q 2 . 

The results for 8a are different in shape and magnitude from those reported by Nagata 



et al. 20|], with the differences more pronounced at large Q 2 . As observed in Figs. H]and[5l 
the dependence on the input form factors and NA couplings is unlikely to account for these 
differences. We have checked the numerical calculations of the TBE amplitudes using two 
independent computer codes, and find agreement between them. It is not clear therefore 
what the origin of the differences may be. Nevertheless, we do agree with the general finding 
in Ref. 20| that the A plays an increasingly important role at forward angles compared with 
the nucleon. 

While the A correction is relatively small for Q 2 between around 0.01 and 3 GeV 2 , 
at very low Q 2 there can be a sizable enhancement of the 7Z contribution at extremely 
forward angles, e — > 1, corresponding to large incident electron energies. This point was 



made recently in Ref. [33J, which argued for a large inelastic Regge contribution in the 
high energy limit. In this region the TPE contribution is suppressed, and the Born term is 
dominated by the proton weak charge, Q w . Hence the A contribution would be enhanced by 



18 



0.04 



0.02 - 



CO 



-0.02 - 



-0.04 - 





= 


■ e 2 


- G 2 


= 0.03 GeV 2 




FIG. 8: TBE correction ^ 7 ( 7 z) arising from nucleon and A intermediate states as a function of the 
incident electron energy E, for Q 2 = (solid) and Q 2 = 0.03 GeV 2 (dashed). 

a factor (1 + Q W )/Q W ~ 14. In Fig. [8] we show the sum of the nucleon and A contributions to 
5 7 ( 7 z) as a function of the incident electron energy, for Q 2 = and for the Qweak 11] value 



J 7 ( 7 Z) 

Q 2 = 0.03 GeV 2 . The A contribution rises linearly with energy up to E ~ 0.5 GeV, where 
it reaches ~ 2 — 3%, after which it decreases. This is qualitatively similar to the resonance 
contributions found in Ref. 



The corrections to the Apy as ym metry at kinematics corresponding to past and planned 
experiments , 111! are listed in Table HJ where the nucleon (8 N ) and A 

(5a) contributions, together with their sum, are shown (in percent %) for various Q 2 and 
laboratory scattering angles 9. In the numerical calculations the empirical proton 25j and 



neutron [27| electromagnetic form factors are used, with dipole parameterizations for the 
axial form factors, as discussed in Sec. IIII1 

For the forward angle HAPPEX Q and GO Q measurements, the nne.eon correction S N 
is in the vicinity of ~ 0.1 — 0.2%, but increases to ~ 1.0 — 1.5% for the backward angle GO 
{4] and the earlier SAMPLE Q| measurements. In contrast, at forward kinematics the A 
contribution 5& is negative and of order —0.5% to —1%, but is almost negligible (~ —0.1%) 
at backward angles. 

When combined, the results reveal a nontrivial interplay between the total nucleon and 
A contributions, with the nucleon dominating the backward angle corrections, and the A 



TABLE I: TBE corrections for the nucleon (<5jv) and A (5a) intermediate states, and their sum 
(in percent), at various experimental kinematics. Also shown are corrections from removing the 



existing hadronic <5^s and total (hadronic + asymptotic) (J™ corrections at Q 2 



Q, 



Q 2 (GeV 2 ) 6 



Expt. 



0.099 
0.477 



6.0° 
12.3° 



HAPPEX [1] 
HAPPEX [1] 



0.077 6.0° HAPPEX [3j 



0.1 144.0° SAMPLE [5] 



0.108 


35.37° 


PVA4 [7] 


1.05 


0.78 


1.83 


0.37 


1.98 


0.23 


35.31° 


PVA4 [7] 


0.62 


0.34 


0.96 


0.23 


1.22 


0.224 


145.0° 


PVA4 [8] 


1.33 


-0.07 


1.27 


0.06 


0.30 


0.122 


6.68° 


GO 


[2] 





18 


-1.06 


-0.88 


0.40 


2.13 
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[2] 
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-0.85 
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_ 
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-0.96 
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0.36 
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-0.92 
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1.85 
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-0.88 
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0.33 


1.77 


0.177 
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[2] 
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-0.83 


-0.67 


0.32 


1.69 
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[2J 
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-0.79 


-0.63 


0.30 


1.60 
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0.95 


0.07 


0.35 
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8.0° 


Qweak [11] 


0.57 
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0.80 


4.25 



5n 5a 5n+a 



0.19 
0.13 



-1.20 -1.01 
-0.44 -0.31 



0.22 



-1.04 -0.82 



1.63 -0.09 1.54 



rhad Atot 
°MS °MS 



0.45 2.42 
0.16 0.86 



0.52 2.78 



0.06 0.33 



contribution driving the forward angle kinematics, where it is rapidly varying with both e 
and Q 2 . Consequently, at the intermediate angles 9 m 35° of the PVA4 experiment both the 
N and A corrections are positive, and combine to give a net ~ 1 — 2% effect. For the planned 



Qweak experiment 



ll| at very low Q 2 (= 0.03 GeV 2 ) and 9 = 8°, on the other hand, the 



positive nucleon and negative A contributions mostly cancel, leaving a much smaller overall 
correction of ~ 0.1%. 

Before correcting the experimental asymmetries for the above TBE effects, one should 



note that the standard data analyses do already include an estimate of 



These are usually taken from the classic analysis of Marciano & Sirlin 
the j(Zj) contributions at Q 2 = 0. Recent explicit calculations 



13 



18 



BE effects [1J,12J]. 



141 ] who computed 
19j . however, have 

found a strong Q 2 dependence at small values of Q 2 , which could significantly impact the 
extrapolation of the Q 2 = results to the experimental kinematics. In order to implement 
the full Q 2 dependence of the TBE corrections, and avoid double counting of the effects 
in the data analyses, one must remove the Q 2 = TBE corrections, which are usually 

nn 

parameterized in terms of p and re pjj, |24J], before adding the corrections computed here. 
In Ref. the loop integration in the box diagram is broken up into a "hadronic", 



low-mass part and an "asymptotic" , high-mass contribution given by 



Ml 



k 2 (k 2 + M 2 ) 



log 



Ml 
p 2 



Mf 



(42) 



where p is the cut-off mass which defines the mass separation, typically of the order of 
1 GeV. For p pa 0.5 — 1 GeV, K asy is in the range ^8 — 10. The hadronic part is computed 
in Ref. [ijj] at Q 2 = using dipole form factors. 

To assess the effect of the new TBE contribution, we display in Table [I] the corrections 
5ms (in percent) defined as 

A v (p, k) - A v (p - Apms, « - Akms) 



S 



MS 



A v (p, re) 

where the numerical values for the ApMS and Akms corrections (for p = 1 GeV) are 



had A _had\ 
^ K MS ) 



( A PMS > 

( a Pms' Ak ms) 



-0.07%, -0.10%) , 
-0.37%, -0.53%) , 



(43) 



(44a) 
(44b) 



for the hadronic only and total (hadronic + asymptotic) 
latter were subtracted in the analyses of Refs. 



18 



contributions, respectively. The 
20j, whereas we believe that only the 
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Q 2 = hadronic component should be removed when adding the new TBE corrections. 
Numerically the hadronic contribution is much smaller than the asymptotic, with the total 
^ms being around 1 — 3% for forward kinematics, and over 4% for the proposed Qweak ex- 



periment 



ll| . The hadronic correction 5^ is also largest at forward angles, but is typically 



0.1 — 0.4% for most of the experiments, and ranging up to 0.8% for the Qweak kinematics. 

The impact of these differences on the strange form factors is difficult to gauge without 
performing a full reanalysis of the data, since in genera 
and form factors are used in the various experiments O, 

an estimate of the induced difference between the strange asymmetry extracted using 
the different form factors was made in Ref. [ljj]. Differences of the order of 15% were found 
between the empirical and monopole form factors (as used in Ref. [3]) for the HAPPEX 



different electroweak parameters 
itE, 7, 8]. Following Zhou et al. 



kinematics 
kinematics 



3|, around 20% for the GO datum in Table I, and over 30% for the PVA4 
7j. One should caution, however, that these values are indicative only, and a 
more detailed reanalysis of the strange form factor data including TBE effects is currently 
in progress [34]. 



V. CONCLUSION 



In this paper we have presented a comprehensive analysis of two-boson (7 and Z) exchange 
corrections in parity- violating electron-proton elastic scattering, paying particular attention 
to the effects arising from the substructure of the nucleon. Working within a hadronic 
framework, we have computed contributions from box (and crossed box) diagrams in which 
the intermediate states are described by nucleons and A baryons. 

The A contribution is found to be much smaller than the nucleon at small e, but becomes 
dominant at forward scattering angles. The dependence of the corrections on the input 
hadronic form factors is small for Q 2 < 1 GeV 2 , but becomes appreciable at higher Q 2 
{Q 2 ~ 5 GeV 2 ), indicating the approximate limit beyond which the hadronic calculations 
may no longer be reliable. 

As well as studying their detailed e and Q 2 dependence, we have evaluated the nucleon 
and A TBE corrections relevant for recent and planned parity-violating experiments 
S y, y, 2, Si, H], finding a nontrivial interplay between the N and A contributions. The 
total corrections at low Q 2 range from ~ —1% for forward angles to ~ 1 — 2% at backward 



22 



kinematics. For the planned Qweak experiment llj we find a large cancellation between 
the (positive) 5n and (negative) 5a corrections, resulting in a modest, ~ 0.1% effect overall. 

Our results for the A differ significantly from those in the recent analysis of Ref. j^oj], 
with the correction 5& differing both in sign and magnitude. We have explored the possible 
origin of these differences by studying the dependence of the corrections on the input nucleon 
and NA transition form factors, but find the effects to be much smaller than that needed to 
explain the discrepancy. We also highlight the need for a careful treatment of the subtraction 
of the standard Marciano-Sirlin 7Z correction at Q 2 = before adding the new contributions. 
The results computed here can be used in future data analyses to more reliably extract 
strange electromagnetic form factors [3] or standard model electroweak parameters 12]. 
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APPENDIX A: RELATIONS TO OTHER iVA TRANSITION FORM FACTORS 

In the literature other notations exist for the NA transition form factors. In this appendix 
we relate the form factors defined in this analysis with those used elsewhere. 



In Ref. 



351 ] (see also Refs. 3l|, |36|) the electromagnetic 7iVA vertex is defined as 

3(M + M A ) [2 f (0 2 )£ ^ 

2M[(M + M A ) 2 + Q 2 ] V 3 \ 9mW ' 

-gEiQ^ip^-g^P-q] ij 5 
gc(Q 2 " 



1 7A— >AT 



(p,q) 



PvQp 



M 



[q 2 (pV 



g 



<f (qy - i a p ■ q)] H5 



To relate this form to that in Eq. (!33j) . we note for the gM term the identity 

e^^upip) = a^ u u a (p) - a m u v (p) + a v V(p) , 



(Al) 



(A2) 



where a^ u = | [7^, 7"], and u a (p) is the Rarita-Schwinger spinor- vector for the spin-3/2 A 
field. Contracting with p v and q a and making use of the constraint relations 

PX(P) = » 7X(P) = > (A3) 
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one finds that the couplings are related by 



9m,e 



9c 



-3M| 



M(M + M A/ 

-3M'i 
M 2 (M + M A ) 



9m,e 



9c 



(A4a) 
(A4b) 



For the axial current, a vertex that one often encounters in the literature is 
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38|) 



32| (see also 



1 1 ZN- 



r A (o 2 \ 



cm 



M 2 



{g a ^p ■ q - q a p^ 



(A5) 



for an outgoing A with momentum p and an incoming Z boson with momentum q. Com- 
paring with the expression in Eq. (1391) . and using the Dirac equation, one finds the following 
relations for the form factors: 



M 
M 2 



2M A 
1 

2M\ 



9?, 



(g£ - 2g£) 



M 2 



q —Q A 

2 93 i 



2M 



(A6a) 
(A6b) 
(A6c) 
(A6d) 

The form factors Cjf and C$ are related by PCAC, C£ —> C A M 2 /Q 2 in the chiral limit, 
with C§(0) = f n gnNA./V3 = 1.2. The fit in Ref. 32| to the neutrino A-production data 
gives C% = and = — leaving a single unique form factor, which is taken to be 
Cg . One may therefore identify the axial couplings in Eq. fl39|) as 



2M 



2 



9t{Q 2 ) 
g£(Q 2 ) 

9i(Q 2 ) 



M 2 



(A7a) 
(A7b) 

(A7c) 



To compute the g£ contribution, we include the 1/q 2 factor in the form factor, and use the 
relation 



q 2 q 2 - A 2 A 2 \ q 2 q 2 - A 
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